Large rapidity gaps survival probabilities at LHC * 
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We calculate the probability of large rapidity gaps in high energy hadronic 
collisions using a model based on QCD mini-jets and soft gluon emission 
down into the infrared region. Comparing with other models we find a 
remarkable agreement among most predictions. 
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1. Introduction 

The possibility of the existence of large rapidity gaps in hadronic colli- 
sions was suggested a number of years ago [HE]- In [T], the occurrence of 
rapidity regions deprived of soft hadronic debris was proposed as a means 
to search for Higgs bosons produced by interactions between colour sin- 
glet particles, e.g. W, Z bosons emitted by initial state quarks, but with 
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a warning against the possibility of soft collisions which would populate 
these regions. For an estimate, an expression for the Large Rapidity Gaps 
Survival Probability (LRGSP) was proposed, namely 

.|o|2._ Jd'h A^^ib,s)\S{braH{b,s) 
' ' Jd^h AAB{b,s)aH{b,s) ^ ' 

where |S'(6)p is the probability in impact parameter space that two hadrons 
A,B got through each other without detectable inelastic interactions, A^^{b, s) 
is the 6— distribution for collisions involved in interactions in which only 
low pt particle emission can take place, and a nib, s) is the cross-section for 
producing, say, a Higgs boson, in such hadronically deprived silent configu- 
ration. To use the above equation, one needs to estimate the probability of 
not having inelastic collisions, |5'(6)p = Pno-inei{b) and the distribution for 
low-pt interactions, namely the probability A^^{b, s) to find those collisions 
for which two hadrons A and B at distance b will not undergo large pt col- 
lisions. Introducing for the scattered partons a cut-off ptmin, above which 
the scattering process can be described by perturbative QCD, we need to 
estimate the 6— distribution of all the collisions with pt < ptmin- In order to 
calculate Pno-ineU can approximate the bulk of hadron-hadron collisions 
with a sum of Poisson distributions for k independent collisions distributed 
around an average n = n(b,s), namely 

oo —k —n 

Pno-^nel = 1 " E ^i^, ^} = 1 " E = ^""^'''^ (2) 

k k=l 

which also leads to the following expressions for the total and inelastic cross- 
sections 

crfrHi = J rf'b[l - e-"(^'^)] (3) 

and 

af,f = 2 J d^h[l - e-"^ cos 3f?ex(6, s)] (4) 

where n(b,s)/2 can be identified with the imaginary part of the eikonal 
function x{b,s). Approximating ?R.ex{b,s) = 0, allows a simple way of cal- 
culating atot if one has a model for n(6, s). A possible strategy to calculate 
the LRGSP is then to build a model for the total cross-section and then in- 
sert the relevant 6— distributions in Eq.([T]). The basic quantity to evaluate 
is thus 

nib, s) = n^'^ib, s) + n'^"^'^(6, s) = ^^^(6, s)c7^p(s) + A'^^^'^ib, s)aw(s) 

(5) 
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where we have spht the average number of colhsions between those with 
outgoing partons below (NP) and above (hard) the cut-ofF ptmin- In the next 
section we shah describe our model [3l d] for total cross-section, comparing 
its results to other models. Then, we shall use the 6— distributions from 
our model to evaluate the LRGSP at LHC, again comparing it to different 
model results. 

2. The Eikonal Mini-jet Model for total cross-section 

To build a realistic model for atot, one needs to understand what makes 
the cross-section rise with a slope compatible with the limits imposed by 
the Proissart theorem, namely atot < log^ s. In our model, the rise in atot 
is driven by the rise of low-x (perturbative) gluon-gluon interactions, while 
the saturation imposed by the Proissart bound comes from initial state 
emission of infrared gluons which temper the too fast rise of the minijet 
cross-section. The rise is calculated using perturbative QCD for collisions 
producing partons with pt > ptmin ~ 1 2GeV, using hard parton scatter- 
ing cross-sections, and the experimentally measured and DGLAP evolved 
parton densities (PDP's) in the scattering hadrons. Thus parton densities 
and the elements of perturbative QCD are the only input needed for the 
calculation of a hard- The rate of rise of this cross-section with energy is de- 
termined by Ptmin and the low-x behaviour of the parton densities. As noted 
before, the rise with energy of the cross-section obtained with this is much 
steeper than that consistent with the Proissart bound, but in our model 
this rise is tempered by soft gluon emission. The saturation mechanism 
takes place through the 5-distribution obtained from the Pourier transform 
of the resummed infrared gluon distribution. This distribution is energy 
dependent and given by [5j 

^2j^^g-iK,.bn(K,) = = ABN{b,qmax) (6) 

where the function h{b,qmax) is obtained through summing soft gluons [5] 
and requires integration of soft gluon momenta from zero to Qmax the max- 
imum transverse momentum allowed by kinematics to single soft gluons. 
The saturation of the Proissart bound is due to the increasing acollinearity 
of "hard" partons produced by initial state soft gluon emission. The sin- 
gle soft gluon distribution needed for the calculation of h{h,qmax) requires 
using infra-red kt gluons and different models with a frozen or singular 
otstrong{kt) producc different saturation effects. We have shown [6j that 
the frozen model is inadequate to quench the rise due to minijets, since 
we see that the early rise in proton- antiproton collisions requires minijets 
with Ptmin ~ 1 GeV , but then the total cross-section rises too much. Put 
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differently, but equivalently, raising ptmin to fit higher energy values of the 
cross-section, say at the Tevatron, would require ptmin ~ 2 GeV but miss 
the early rise. Instead, we find that a singular a<j produces an adequate 
s— dependent saturation effect. Singular expressions for Os are discussed in 
the literature, in particular for quarkonium phenomenology [7J. Our choice 
is a singular, but integrable expression for the strong coupling constant in 
the infrared region, namely 

^ ^ 33 - 2Nf^ kt ' ^ ^ 

where the scale factor is chosen to allow a smooth interpolation to the 
asymptotic freedom expression for as, namely we choose 

127r p 
*^ ^ 33-2iV^log(l+p(^)2p) ^ > 

The singularity in the infrared is regulated by the parameter p, which has 
to be < 1 for the integral in h{b,qmax) to converge. The next input for 
phenomenological tests of our model is the number of non-perturbative (NP) 
collisions. We approximate it as 

n^^ = A%a,{l + ^) (9) 

with e = 0, 1 for the process pp or pp. We choose a constant fio ~ 48 mb 
and use for the same model as for the hard collisions, but we restrict 
Qmax to be no larger than ~ 20%ptmm) since these collisions are limited to 
Pt < Ptmin- The same function is then used for the LRGSP calculation, 
as discussed in the next section where both the estimated atotai as well as 
the LRGSP will be presented and compared with other models. 



3. Total cross-sections and survival probability 

Applying the above described model to the calculation of total cross- 
sections, gives the results shown in the left panel of Figure [TJ To obtain 
this figure, we have used different PDF's and slightly different values for the 
parameters p, ao and ptmin and the variations are indicated by the band. We 
have shown [3] that the asymptotic behaviour of this cross-section can be 
fitted with a log^ s type behaviour. This is a phenomenological confirmation 
that the model satisfies Froissart bound. 

We do not see in our model any hard Pomeron behaviour beyond the 
initial rise and predict a value = lOOtllimb). The LRGSP can now 

be calculated using the two quantities e""^^'''^ and A^^{b,s) which were 
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Fig. 1. Total cross-section data and models [9] (left) and survival probability for 
large rapidity gaps for different models (right). 



input in the calculation of the total cross-section. The right panel of Figure 
[T] shows such evaluations, with the yellow band corresponding to various 
MRST densities [TT], the central full line to using GRV densities |12| and 
the other lines and bands representing comparisons with other models, as 
indicated. The most interesting result of this figure is that the predictions 
for LHC agree reasonably well amongst eachother, namely \S\^ = 5 10%, 
in spite of the fact that the various models differ greatly in details and the 
way in which they achieve results for total cross-sections consistent with the 
Froissart bound. Thus the model estimates for the LRGSP are quite robust. 



4. Conclusions 

We have built a model for utot which incorporates hard and soft gluon 
effects, satisfies the Froissart bound and can be used reliably to study other 
minimum bias effects e.g. the Survival Probability of Large Rapidity Gaps. 
It can also be extended to calculations of total 7p and and 77 cross-sections. 
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